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PREFACE 


Tills  report  is  an  account  of  the  work  performed  by  the  McDonnell 
Douglas  Research  Laboratories  on  Laser  Pulse  Speckle  Effects  for  the  U.S. 
Army  Missile  Research  and  Development  Command  (MIRADCOM)  under  subcontract 
to  tlie  Battelle  Columbus  Laboratories  (Scientific  Services  Agreement 
D.O.  No.  0522)  and  authorized  by  the  Army  Research  Office  Contract 
DAAG29-76-D-0100 . The  work  was  performed  in  the  Chemical  Physics  Labora- 
tory, managed  by  Dr.  C.  J,  Wolf.  The  principal  investigator  was  Dr. 

J.  C.  Leader;  Mr.  J.  M.  Putnam  was  responsible  for  the  computer  programming. 
The  contracting  officer's  technical  representative  was  Dr.  John  P.  Stettler, 
Physical  Sciences  Directorate,  MIRADCOM,  Huntsville,  AL. 
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i.  OBJECTIVE  AND  AI’PKOACH 


pmm 


The  objective  of  this  study  was  to  investigate  the  temporal  correlations 
of  laser  light  scattered  by  a target  illuminated  with  a laser  pulse  and 
assess  tlie  effects  of  these  correlations  on  the  rime/speckle-spectrum  signa- 
ture of  the  target.  Because  of  tlie  limited  duration  and  level-of-ef fort 
provided  by  the  Scientific  Services  Agreement,  the  investigative  approach 
focused  upon  analyzing  the  phenomena  from  an  overall  viewpoint  and  utilized 
previous  analytical  and  computational  tools  maximally.  To  achieve  maximum 
benefits  to  MIRADCOM  in  the  allocated  time  period,  emphasis  was  placed  upon 
(a)  a physical  understanding  of  the  scattering  phenomena,  (b)  a qualitative 
assessment  of  the  impact  of  pulse  speckle  effects  upon  pulse-Doppler  LADAR 
systems,  and  (c)  the  identification  of  appropriate  analytical  methods  v/hich 
can  be  used  for  more  detailed  analyses  of  the  statistics  of  time/speckle- 
spectrum  target  images. 
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2.  SUMMARY  OF  RESULTS 


Significant  results  of  this  contract  study  are  sununarized  below,  A 
more  detailed  outline  of  analytical  and  computed  results  is  provided  in  Sec- 
tions A and  5. 

• Target  rotation  during  a pulse  interval  degrades  the  performance  of 
a linear  FM  chirp  LADAR  range-Doppler  system  because  many  speckle 
lobes  are  swept  by  the  receiver. 

• Intra-pulse  speckle  can  degrade  the  performance  of  coherent  pulse- 
train  LADAR  range-Doppler  systems  because  time  samples  of  the  spectra 
are  smeared  in  range  by  speckle  fluctuations.  A priori  knowledge  of 
the  target's  range  and  velocity  can  reduce  this  degradation  via 
suitable  data  processing. 

• The  statistical  fluctuations  of  laser  light  scattered  by  a three- 
dimensional  target  illuminated  with  a laser  pulse  are  generally  non- 
stationary. However,  the  fluctuations  can  be  treated  as  a random 
process  with  stationary  increments  and  described  by  an  appropriately 
defined  structure  function.  Amplitude  statistics  are  generally 
zero-mean,  circular  Gaussian. 

• An  analytical  method  for  calculating  (ensemble-averaged)  undistorted 
range/speckle-spectrum  signatures  and  structure  functions  of  conical 
targets  was  developed.  Calculational  limitations  imposed  by 
customary,  flat-plate,  target  quantization  techniques  are  obviated 
by  this  method. 

• The  surface  roughness  of  conical  targets  impacts  the  LADAR  time/ 
speckle-spectrum  signature  most  dramatically  for  laser-beam  incident 
angles  approaching  the  (cone-axis)  normal. 

• The  scattered-field  fluctuations  resulting  from  a laser-pulse 
illuminated,  homogeneous,  randomly  rough  cone  have  normalized  local 
homogeneity  (as  defined  in  this  investigation). 
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# Detailed,  parametric  caicuiations  of  speckle  distortion  effects  should 
be  performed  to  assess  the  maximum  mensurational  capability  (i.e., 
ignoring  other  noise  effects)  of  both  chirp  and  coherent  pulse-train 
time/speckle-spectrum  LADAR  target  imaging  systems.  Results  of  these 
calculations  should  be  compared  with  analogous  experimental  data,  as 
available,  to  guide  subsequent  analytical  and  experimental  efforts. 

# Because  the  scattered-field  structure  function  is  central  to  calcula- 
tions of  range-smearing  effects  for  coherent  pulse-train  systems,  the 
structure  function  for  non-conical  targets  (of  interest  to  BMD  and 
MIRADCOM)  should  be  performed  to  test  the  similarity  of  these  structure 
functions  with  the  functional  form  (see  Section  5)  found  for  conical 
targets. 


4.  ANALYSIS 


4. I Pulse  Speckle  Effects  on  Tlme/Speckle-Spectrum  Signatures 

Basically,  two  types  of  range-Doppler  imaging  systems  have  been  proposed 
to  mensurate  targets  via  their  time/speckle-spectrum  signature.  The  linear 
FM  chirp  system^  determines  the  target  range  by  demodulating  the  received 
(scattered)  chirp  pulse  so  that  range  Is  proportional  to  the  frequency  of  the 
demodulated  signal.  Doppler  (spc.'kle  spectrum)  information  is  obtained  by 
Fourier  transforming  numerous  pulse  amplitudes  for  each  range  cell.  Receiver 
bandwidth  requirements  are  generally  not  stressed  by  this  technique.  The 
coherent  pulse-train  technique  relies  upon  the  limited  time  duration  of  the 
pulses  comprising  the  pulse-train  to  range-resolve  the  target  while  the  spec- 
trum of  fluctuations  is  sampled  by  the  pulses  so  that  the  pulse-train 
duration  determines  the  spectrum  resolution.  The  receiving  system  bandwidth 
is  generally  determined  by  the  inverse  pulse-width  for  the  coherent  pulse- 
train  system.  The  large  bandwidths  resulting  from  detailed  range  resolution 
requirements  stress  the  current  optical  receiver  technology. 

Random,  radiation  field-fluctuations  (producing  the  irradiance  fluctua- 
tions commonly  known  as  speckle)  result  from  basically  two  types  of  laser 
scattering  phenomena:  (1)  rough-surface  laser  scattering  reduces  the  spatial 

coherence  of  the  laser  light  such  that  the  scattered  field  is  only  partially 

2 

correlated  at  different  scattering  angles,"  and  (2)  the  illumination  of 
different  rough-surface  areas  decorrelates  the  scattered  field  because  the 
rough-surface  areas  are  only  partially  correlated.  Similar  coherence  losses 

3 

result  from  random  volume-scattering  processes.  The  coherence  losses  re- 
sulting from  scattering  angle  differences  and  different  illuminated-areas 
impact,  respectively,  the  linear  FM  chirp  and  coherent  pulse-train  imaging 
systems . 

Because  the  pulse  width  of  the  linear  FM  chirp  range-Doppler  imaging 
system  is  generally  much  longer  than  any  target  dimension  of  interest  (L/c, 
where  L is  the  target  length  and  c is  the  velocity  of  light),  coherence 
losses  resulting  from  different  illuminated  areas  do  not  significantly 
influence  this  system.  However,  target  rotation  during  the  time  the  pulse 
illuminates  the  target  can  produce  a temporal  amplitude  fluctuation  of  the 
scattered  pulse  resulting  from  the  changing  scattering  coordinates  of  the  re- 
ceiver. The  number  of  spatially  de-correlated  (i.e.,  speckle)  cells  seen  by 


the  receiver  during  the  pulse  Interval  (T^^)  is  given  approximately  by 


N = 2QT  L/A  , 
s p 


where  is  the  target  rotation  rate  and  A is  the  laser  wavelength.  The  fre- 

4 

quancy  uncertainty  resulting  from  these  fluctuations  is 


Av  = 2Ln/A  . (2) 

This  frequency  uncertainty  can  be  comparable  with  the  frequency  range  deter- 
li'ining  the  length  of  the  target,  and  thus  the  range-Doppler  image  of  the 
target  is  degraded. 

Because  the  pulse-width  of  the  coherent  pulse-train,  range-Doppler 
imaging  system  is  much  less  than  the  range  dimension  of  the  target  (to  range 
resolve  it),  the  number  of  speckle  cells,  subtended  during  a pulse 
interval  is  fractional  for  most  rotation  rates  of  interest  at  optical  wave- 
lengths. Thus,  coherence  losses  resulting  from  scattering  angle  changes, 
during  a pulse  interval,  are  negligible.  However,  because  the  pulse  interval 
is  less  than  the  target  dimension,  independent  rough-surface  areas  are 
illuminated  as  the  pulse  traverses  the  target,  and  the  resulting  scattered 
field  fluctuations  can  alter  the  expected  range-Doppler  image.  Because 
MIRADCOM  is  primarily  interested  in  the  coherent,  pulse-train,  range-Doppler 
imaging  system,  the  remainder  of  this  report  is  devoted  to  analysis  of  the 
temporal-field  fluctuation  resulting  from  the  laser  illumination  of  uncor- 
related surface  areas  and  their  Impact  upon  the  time/speckle-spectrum 
signature. 

For  simplicity  the  square-wave  pulse  train  illustrated  in  Figure  1 is 
considered,  and  the  classic  BMD  conical  target  (also  illustrated)  is  treated. 
Denoting  the  incident  electric  field  by  E (t)  and  the  field  scattered  by  the 
target  to  point  P (at  a range,  R,  from  the  target)  with  polarization  compo- 
nent j by  E ,(t),  the  average  scattered -f ield  intensity  is  (by  definition) 
s J 

o(t-R/c), 


where  0(t;)  is  the  cross  section  of  the  target  that  is  illuminated  by  the 
pulse  at  time,  t.  Speckle  effects  resulting  from  uncorrelated,  illuminated 
target  areas  cause  the  actual  intensity  to  fluctuate  about  the  mean  value 
predicted  by  Equation  (3).  Thus,  for  a given  pulse  comprising  the  pulse 
train,  tlie  magnitude  of  the  actual  scattered  field  could  fluctuate  in  a 
manner  similar  to  the  fluctuations  shown  in  Figure  2.  Because  the  variance 
of  che  field  fluctuations  is  dependent  upon  time,  the  statistics  of  the 
field  fluctuations  are  clearly  non-stationary . Arguments  similar  to  those 
presented  in  Reference  4 indicate  that  the  field  becomes  decorrelated  in  a 
time  interval  on  the  order  of  a pulse  width,  T . 
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Figure  1 Pulie-train  parameters 
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If  the  scattered  field  Is  coherently  detected,  the  Intermediate 
frequency  (IF)  photocurrent  output,  x(t),  is  proportional  to  the  scattered 
field  (neglecting  finite  aperture  effects,  shot-noise,  and  finite  bandwidth 
filtering  effects).  Thus,  assuming  that  the  range  rate  of  the  target  is 
zero,  the  ideal  (i.e.,  uncorrupted  by  noise  and  filtering  effects)  photo- 
current signal  varies  temporally  in  a manner  similar  to  that  shown  in  Figure 
3.  The  speckle  spectrum  within  a range  interval  corresponding  to  a pulse- 
width  (Tp)  is  obtained  by  squaring  the  Fourier  transform  of  some  measure  of 
the  field  (photocurrent)  fluctuations  resulting  from  scattering  within  that 
range  interval.  A reasonable  measure  of  the  field  is  provided  by  the 
integrated  photocurrent  output  within  a pulse  interval,  i.e.. 


X(t  ,AR) 
n 


t = 2(R  + AR) /c  + (n  - 1)T^ 


= 2(R  + AR) /c  + (n  - 1)T^ 


+ T 

P 


(4) 


t 


ti=2l 


t 


t2  “ 2R/c  + 
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The  sampled  (discrete)  Fourier  transform  of  the  field  fluctuations  is  then 


(u.AR) 

c 


Si 

E 


X(t^,AR)e"^“*^n 


for  the  N pulses  comprising  the  pulse-train  and  the  square-law  spectrum 
(proportional  to  the  target's  cross-range,  differential  cross  section)  is 


S(W.AR)  = Ix^  (W.AR)!^, 

c c 


where  is  the  coherence  time  (l.e.,  duration  of  the  pulse-train).  If  the 
nearest  part  of  the  target  is  initially  at  a range  R (at  the  beginning  of  the 
pulse-train)  and  has  a relative  velocity. 


R - ^ 

^ dt  ‘ 


the  appropriate  measure  of  the  field  fluctuations  resulting  from  the  range 
increment  at  a range  AR  along  the  target  [corresponding  to  the  static  measure 
provided  by  Equation  (4) J is 


X<t„.r) 


t - 2(E  + AR)/c  + (n  - 1)1  (1  + R/c)  + T 

f ' ' 

I dtx( t) 

t = 2(R  + AR)/c  + (n  - 1)T^(1  + R/c) 


The  correction  to  the  integration  Interval  [(n  - l)T^R/cJ  shown  in  Equation 

(8)  is  negligible  for  most  relative  velocities  and  coherent  integration  times 

[T  = (n  - 1)T  1 of  Interest  because  of  the  small  value  of  the  ratio  R/c. 
c r 

However,  the  cumulative  effect  of  this  time  (range)  interval  correction 
(T^R/c)  over  the  Incoherent  Integration  time  (T^) , required  to  reduce  the 
variance  of  the  deduced  speckle  (Doppler)  spectrum,  is  significant  because 


T >>T  . 
i c 
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An  examination  of  the  systems  consequences  of  intra-pulse  speckle  does 
not  fall  within  the  purview  of  this  investigation.  However,  because  the 
time/speckle-spectrum  signature  consequences  of  Intru-pulse  speckle  are 
system  dependent,  it  is  worthwhile  to  itemize  a few  examples  of  the  types  of 
effects  that  can  be  expected. 

(1)  If  the  receiving  system  is  designed  to  time  (range)  the  target  by 

sensing  some  measure  of  the  arrival  time  of  each  pulse  (e.g.,  leading  edge, 

trailing  edge,  mean,  etc.),  inaccuracies  in  the  arrival  time  on  the  order  of 

a pulse  width  will  result  from  speckle  fluctuations.  Thus,  a ^andom  mixture 

of  adjacent  range  increment  measures  of  the  photocurrent  [X(t^,AR  + R^) 

where  R is  a random  variable  with  a variance  on  the  order  of  cT  ] are 
n p 

Fourier  transformed  to  obtain  the  cross-range  spectrum,  S(ll),AR).  The  re- 
sultant time/speckle-spectrum  signature  is  consequently  degraded. 

(2)  If  the  receiving  system  establishes  fixed  range-interval  increments 
to  sample  field  fluctuations  [in  the  manner  suggested  by  Equation  (4)], 
admixtures  of  adjacent  range  intervals  [as  suggested  in  (1)]  in  the 
discrete  Fourier  transform  of  Equation  (5)  can  result  from  receiving  system 
jitter. 

(3)  If  t.ie  target's  relative  velocity  is  unknown  (or  improperly 
estimated)  and  fixed  range  Increments  are  used  [as  suggested  in  (2)],  the 
incoherent  integration  process  will  include  square-law  spectrum  samples  from 
adjacent  range  intervals.  Although  the  sample  spectra  are  quickly  de- 
correlated  in  range,  the  resultant  averaged  spectrum  has  a degraded  range 
and  frequency  resolution. 

The  examples  cited  above  illustrate  the  manner  in  which  intra-pulse 
speckle  can  influence  time/speckle-spectrum  signatures.  The  precise  manner 
in  which  speckle  affects  the  signature  derived  from  a particular  range- 
Doppler  imaging  system  is  obviously  dependent  upon  the  details  of  the  system 
implementation.  To  analyze  speckle  effects  on  any  system,  however,  it  is 
necessary  to  specify  the  statistics  of  the  speckle  fluctuations.  The 
statistics  of  intra-pulse  speckle  are  explored  below. 


4 • 2 Pulse  Speckle  Statistics 

The  probability  distribution  function  (PDF)  describing  speckle  fluctua- 
tions can  be  quite  generally  described  by  the  Hoyt  distribution.^  However, 
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because  the  coherent  ccntrlbution  to  the  scattered  field  results  from 
approximately  equidistant  scattering  centers  to  achieve  a coherent  super- 
position (i  e. , glint  returns),  the  number  of  target  glints  is  usually  small 
and  spatially  limited.  Therefore,  the  ensemble  average  PDF  of  the  field 
amplitude  for  the  majority  of  the  time  domain  of  the  scattered  pi.l8e  must 

g 

have  non-circular  Gaussian  statistics.  In  fact,  because  the  scattered 
signal  results  from  rough-surface  scattering  at  various  ranges  considerably 
in  excess  of  a laser  wavelength,  the  PDF  for  each  polarization  component  re- 
duces to  circular  Gaussian  statistics. 

Because  the  statistics  of  the  scattered  pulse  amplitude  are  circular 
Gaussian  (i.e.,  the  real  and  imaginary  parts  have  equal  variances),  it  is 
sufficient  to  describe  the  correlation  function  of  the  field  to  completely 

9 

specify  the  statistical  description  of  field  fluctuations.  However,  as 
noted  earlier,  the  temporal  field  amplitude  fluctuations  are  statistically 
non-stationary . Therefore,  the  correlation  function  is  an  inappropriate 
statistical  descriptor  of  the  field  amplitude  fluctuations.  However,  because 
(a)  the  scattered  field  amplitude  increment, (t)  - E^j(t')^  Is  time  in- 
dependent (t)^  “ ^E^j  (t  ' )^  * oj  and  (b)  a suitably  defined  structure 

function  (shown  below)  depends  only  on  the  time  differences  t - t/2  and 
t + t/2,  the  field  fluctuations  can  be  treated  as  a random  process  with 
stationary  increments . The  appropriate  structure  function  for  studying 
fluctuations  of  any  combination  of  polarization  components  of  the  scattered 
field  is 


JDg^'Ct.T)  + T/2)  - E^^(t  - T/2)][Egj^(t  + T/2)  - E^j^Ct  - t/2)1*) 

5^E^j(t  + T/2)E*^(t  + t/2)^  +<^E^j(t  - T/2)E*^(t  - t/2))>  (9) 

-<Esj(t  + T/2)E*^(t  - T/2>-<E^j(t  - T/2)E*^(t  -H  t/2)>. 

The  structure  function  defined  by  Equation  (9)  characterizes  the  intensity  of 

fluctuations  having  a period  f t.  The  structure  function  can  also  be  ex- 
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pressed  in  terms  of  the  power  spectrum  of  field  fluctuations,  (u) , as 


mmmiiswii^: 


j!lj>.ijPi,Wi..«>il»^>iJilM.UU-  IJA-^  W,  Jlgy 


j 


cos  (i)(t2 


t,)]Jw^*^(co)d(o. 


(10) 


Thus,  if  the  power  spectrum  of  field  fluctuations  is  charactet ized  by  e 
2 

sine  (u))  function,  l.e., 


(11) 


the  structure  function  is  given  by 


, ( 7r/2  t/2T  t ^ 2T 

\ = { . (12) 
I Tr/2  T > 2T 


Because  the  scattered  field  consists  largely  of  the  incoherent  component  (by 
the  above  arguments),  correlations  at  different  time  Increments  result  from 
overlapping  illumination  areas  on  the  target  (see  Figure  1).  Therefore,  the. 
Equation  (9)  definition  reduces  to 


^D^^(t,T)  =<(Egj(t  + 'r/2)E*^(t  + T/2)^+<^Egj(t  - T/2)E*j^(t  - t/2)^ 
- 2<^gj(t  + T/2)E*j^(t  - T/2)^. 


(13) 


Since  the  scattered  field  intensities  comprising  the  right-havid  side  of 
Equation  (13)  are  proportional  to  an  appropriate  range-resolved  target  cross 
section  [see  Equation  (3)],  the  structure  function  defined  by  Equation  (13) 
is  clearly  a function  of  time,  t.  The  scattered  field  is  therefore  a non- 
homogeneous  random  process  with  stationary  Increments.  However,  if  an 
appropriately  normalized  structure  is  defined,  viz.. 


(lA) 


and  if  the  normalized  structure  function  is  largely  time  independent  within 
the  scattered  pulse  duration,  i.e.. 
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(15) 


^ for  2R/c  - T /2  < t < 2(R  + L)/.'.  + T /2  , 

6 © P r 

the  utility  of  the  structure  function  analysis  for  describing  pulse  amplitude 
fluctuations  is  considerably  enhanced.  Field  fluctuations  having  structure 
functions  that  satisfy  the  Equation  (15)  relation  are  defined  in  this  inves- 
tigation as  having  normalized  local  homogeneity  in  analogy  with  the  standard 
definition  of  locally  homogeneous  random  processes.  The  increased  utility  of 
the  analysis  results  from  the  observation  that  for  most  targets  of  interest, 

+ T/2)) 

- T/2)E*^(t  - 1/2))]  (1,) 

for  T < 

by  the  mean  value  theorem,  providing  the  cross  section  is  sufficiently 

jmoothly  varying  and  T is  not  too  large.  Thus,  within  the  range  of 

s 

validity  of  the  Equation  (16)  approximation,  a function  analogous  to  the 
degree  of  coherence  can  be  defined,  i.e.. 


-k  T/2)E;^(t  - T/2)) 


(17) 


and  expressed  in  terms  of  the  normalized  structure  function,  i.e., 


Yj°(t,T)  1 - ^a)^(t,T) 


(18) 


for  T < Tg. 


The  function  defined  by  Equation  (17)  describes  the  local  correlations  of  the 
scattered  field  and  has  the  normalized  properties 


0 - Yjj°(t,T)  ::  1 

YjJ(t.O)  E 1 


(19) 


'f>r^ 
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akin  to  the  conventional  degree  of  coherence.  Therefore,  this  function  will 
be  called  the  local  degree  of  coherence.  Clearly,  if  the  scattered  field  lias 
normalized  local  homogeneity  [Equation  (15)  is  valid],  the  local  degree  of 
coherence  Is  temporally  invariant,  l.e., 

YjJ(t,T)  ^ YjJ(t)  for  2R/c  - 1^/?.  < t < 2(R  + L)/c  + 1^/2.  (20) 

This  Important  property  permits  the  calculation  of  the  deleterious  intra- 
pulse speckle  effects  (discussed  in  Section  4.1)  from  a knowledge  of  the 
uncorrupted  time/speckle-spectrum  signature  and  the  normalized  structure 
function.  The  remainder  of  this  section  examines  means  for  calculating  the 
structure  function,  normalized  structure  function,  and  (as  an  interesting 
byproduct)  the  various  resolved  signatures  of  a conical  target  having 
various  surface-roughness  properties.  Section  5 provides  calculations  that 
demonstrate  that  conical  targets  indeed  have  normalized  local  homogeneity. 

4. 3 Conical  Target  signature  and  Structure  Function  Calculations 

If  a pUr  e-wave  laser  pulse-train  illuminates  a conical  target,  the 
intersection  of  the  planar  wave-front  with  the  cone  describes  a conic  section 
(i.e.,  truncated  ellipse  or  hyperbola).  The  time/speckle-spectrum  signature 
is  proportional  to  the  differential  laser  radar  cross  section  arising  from 
cross-range  elements  within  the  time  period  of  the  pulse.  The  differential 
cross  section  of  the  cross-hatched  scattering  area  illustrated  in  Figure  4 
is  proportional  to  the  time/speckle-spectrum  signature  of  the  conical  target 
for  a time  delay  2t  and  speckle  spectrum  at  x having  a differential  spectral 
width  dx,  for  a plane  wave  incident  at  the  polar  angle  0 (in  the  y-z  plane). 
Normalized  spectral  units  of  x and  do/dx  are  used  in  this  investigation  as 
in  Reference  12  so  that  knowledge  of  the  cone  radius,  spin  rate,  range, 
incident  angle,  wavelength,  and  laser  power  permits  the  calculation  of  the 
absolute  frequency  spectrum  in  units  of  power  and  frequency.  The  projection 
of  the  plane-wave  conic  intersection  onto  tlie  base  of  the  cone  also  ^escribes 
a conic  section.  The  base  plane  (i.e.,  x-y)  coordinates  of  the  time  and 
spectrum  intersections  can  be  solved  using  the  equation  of  the  plane 


I 

I 

I 


! » 

! -i 


z = tan  0 y -1  ^ , 


(21) 


14 


The  intersection  of  the  cone  and  plane  is  found  from  equating  Equations 
(21)  and  (22),  i.e., 


-tan  0(Y  - a)  ± ^ - a)^  + a^x  ^ tan^Q  - a^x  ^ 

' n n 

(tan^Q  - a^) 


(23) 


where  Y - i/R^  and  x^  - x/R^.,  » y/Rt  are  normalized  units  having  maximum 

magnitude  of  unity.  The  requirement  that  the  Intersection  lies  on  the 
finite  cone  rather  than  the  infinite  cone  is  mathematically  stated 


< 


1. 


(24) 


The  differential  area,  A^.,  of  the  spectrum-pulse  width  intersection  can  be 
found  from  the  relation 


Aj.  Ag  csec  3, 


(25) 


where  3 is  th'  cone  half-angle  and  Ag  is  the  base  projected  area.  The  base 
projected  area  can  be  solved  analytically  from  the  expression 


+ c dx 


~~—y 

+ c 

1 

+ c 


+ -■  - log  (xV7  +Vax^  + c')  for  a > 0 
2Va 

(26) 

+ — - — sin  ^(xV^)  for  a < Q 

2V^ 


where  a and  c are  determined  from  the  major  and  minor  axes  of  the  elipse  or 
hyperbola  projected  upon  the  base.  The  differential  cross  section  of  the 
pulse-spectrum  intersection  is  then 


da(t,x,Tp) 

dx 

P 


- Aa(t,x,T^) I 

P.  prnj 

Ax 

P 


(27) 


where  a^(a))  is  the  laser  cross-section  per  unit-area  of  the  differential 
element  at  the  azimuthal  angle  oj.  The  average  azimuthal  angle,  id,  can  be 
found  from  the  relation 


CD(x) 


•j[a)(x,t) 


+ u)(x,  t + 


T )]. 

P 


(28) 


where 


cj(x,t)  - tan  ^fx  /y  (x  , t)].  (29) 

n n n 

The  laser  cross-section  per  unit-area,  C°,  can  be  calculated  from  the  com- 
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puter  code  ROSSCO,  using  the  relations  in  Reference  12  relating  the  local 
cone  incident  and  polarization  angles  to  the  cone  laser  incident  angle  and 
polarization  angles  via  the  azimuthal  angle.  Cone  shadowing  effects  are 
accounted  for  via  the  caveats 


da(t,x,T  ) 1 

E_lproj  = Q 

dx 


(30) 


if 


u(x,t)  > U)^  and  a)(x,t  + T^)  > 


(31) 


and 
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i| 

■t  i 


da(t,x,T 


dx 


P pro1 


AjXt,x,Tp)0°(w)  F(w), 


(3 


if 


and 


t*)(x,t) 


< 10^ 


(3. 


a)(x,t  + T ) > ui, , 


where 


F(u)) 


- a)(x,t)  |^|ai(x,t)  - u)(x,t  + T )j, 


(3^) 


and  0^  is  the  shadowing  azimuthal  angle  defined  in  Reference  12. 

Equations  (23)  - (34)  and  the  requisite  Equations  of  Reference  12  when 
computer  implemented  with  suitable  logic  required  to  handle  special  integra- 
tion cases  arising  near  the  edge  of  the  cone,  permit  accurate  calculations 
of  the  ideal  time/speckle-spectrum  signature  for  arbitrary  scattering 
geometries,  pulse  widths  and  cone  roughness  conditions.  Significantly,  this 
analytical  method  of  calculating  signatures  does  not  require  subdividing  the 
target  into  approximate  flat-plate  areas.  Therefore,  small  pulse-width 
signatures  (required  for  subsequent  structure-function  calculations)  are 
facilitated,  computer  storage  requirements  are  minimal,  and  the  required 
calculations  are  extremely  rapid.  The  time  (range)  resolved  cross  section  is 
given  by 


o(t,T  ) , 

' p' 'proj 


da(t,x,T  ) . 

dx  • * * p " ' proj . 

dx 


(35) 


Therefore  the  structure  function  for  the  cone  is  given  by 


^D^(t,T)  = K[o(t  + t/2,T  ;j,k)  + a(t  - T/2,T^;j,k) 
e P P 


- 2o(t,Tp  - r;j  ,k)]  , 


(36) 


where 


4itR 


and  the  normalized  structure  function  is 
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(37) 


^!jOg(t,T)  = [o(t  + T/2,Tp;j,k)  + a(t  - T/2,Tp;j,k) 

- 2a(t,Tp  - T;j ,k)]/o(t,Tp;j ,k) , 

where  the  argument  notation  j and  k has  been  added  to  denote  the  polarization 
components. 

Section  5 provides  sample  calculations  of  tlme/speckle-spectrum  signa- 
tures [Equation  (27)],  time-resolved  cross-section  distributions  [Equation 
(35)],  structure  functions  [Equation  (36)]  and  normalized  structure  functions 
[Equation  (37)]  for  various  cone  incident  angles  and  target  roughness'^s . 
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5 . COMPUTED  RESULTS 


The  expressions  outlined  in  Section  4.3  were  computer  coded  to  obtain 
numerical  results  for  the  signatures  and  structure  functions,  discussed  in 
Sections  4.1  and  4.2,  for  a conical  target.  A cone  with  a 10°  half  angle  was 
assumed,  and  a normalized  base  radius  of  unity  was  used  in  accord  with  the 
normalization  (scaling)  procedure  discussed  in  Reference  12.  All  calcula- 
Lions  assumed  incident  and  scattered  polarization  vectors  that  were  normal  to 
the  plane  of  incidence  to  the  cone  axis  (i.e.,  HH) . Four  laser  Incident 
angles  with  respect  to  the  cone  axis  (0  ■ O',  45°,  60°,  89°)  and  two  rough- 
ness conditions  were  examined.  A nominal  surface-roughness  condition 

corresponding  to  phenolic  carbon  scattering  at  10.6  pm  was  employed  which  was 
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determined  from  a ROSSCO  analysis  of  actual  laser  scattering  data.  The 

derived  ROSSCO  parameters  are 

Normal  Phenolic  Carbon 

ko  ■=■  0.529,  ki  - 2.72,  s - 0.246,  m « 4.0,  n * 1.49,  a - 6.47  x 10“^. 

A smoother  phenolic  carbon  material  was  examined  with  the  following  ROSSCO 
parameters: 

Smooth  Phenolic  Carbon 

ko  ■ 0.529,  ki  ■ 5.0,  s * 0.1,  m ■ 4.0,  n * 1.49,  a “ 6.47  x 10 
A pulse  width,  T^,  of  0.33  (in  units  of  the  normalized  radius,  x)  was  assumed, 
and  normalized  time  intervals  beginning  at  zero  at  the  nearest  range  of  the 
conical  target  were  used.  Computed  results  are  presented  and  discussed  in 
the  following  subsections. 

5 . 1 Time/Speckle-Spectrum  Signatures 

Computed,  ideal  (i.e.,  noise-free)  time/speckle-spectrum  signatures  are 
illustrated  three-dlmensionally  in  Figures  5,  6,  7,  and  8 for  the  0°,  45°, 

60°,  and  89°  incident  angles,  respectively.  Only  positive  values  of  x are 
depicted  because  the  calculated  signature  is  symmetric  (i.e.,  the  values  of 
dc/dx  for  negative  x values  are  mirror  images  of  the  values  for  positive  x) . 
The  spectral  width  illustrated  in  Figure  5 for  the  0°  incident  angle  is 
somewhat  fictitious  because  a normalization  constant  of  zero  (sin  0) 
multiplies  the  x axis.  Nevertheless,  a signature  similar  to  the  computed  one 
is  expected  for  nearby  incident  angles. 
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Two  general  observations  concerning  the  calculated  signatures  can  be 
made.  First,  It  Is  noteworthy  that  for  Incident  angles  near  nose-on  (G  - 0®), 
the  increased  scattering  area  from  the  edge  of  the  cone  enhances  the  edge 
definition.  This  effect  persists  to  an  incident  angle  of  60®  for  carbon 
phenolic  of  normal  roughness.  Secondly,  the  effects  of  surface  roughness 
(other  than  affecting  the  cross-section  magnitude)  become  more  pronounced  for 
larger  angles  of  incidence  with  respect  to  the  cone  axis.  The  aforementioned 
edge  enhancement  (resulting  from  the  increased  scattering  area)  is 
effectively  eliminated  by  the  smoother  phenolic  carbon  for  an  incident  angle 
of  45®.  The  numerical  results  indicate  that  it  is  doubtful  if  the  edge  of 
the  smooth  target  can  be  resolved  for  angles  of  incidence  of  60®  or  greater. 

5.2  Time-Resolved  Cross  Sections 

Computations  of  the  time-resolved  cross  section  [Equation  (35)]  for  the 
four  incident  angles  of  interest  are  illustrated  in  Figures  9-12  for  a 
time  resolution  of  At  - 0.05  (normalized  units).  It  is  noteworthy  that 
deviations  from  a ramp  distribution  are  observed  as  the  incident  angle  in- 
creases and  that  the  surface  roughness  alters  the  shape  of  the  distribution 
for  the  larger  incident  angles. 


0-0® 


Time  (normalized  units) 


Figure  0 Time-resolved  cross  section  for  phenolic  carbon  cone  10  - 0®) 

aerr-sttr-r 
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Time  (normalizud  units) 

Figure  10  Time-resolved  cross  sections  for  phenolic  carbon  cone  45°) 
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5.3  Cone  Structure  Functions 

2 

Calculated  cone-structure  functions  [Equation  (34)  In  units  of  LRCS/R^ 
and  T/Tp]  are  Illustrated  In  Figures  13  - 16  for  the  four  Incident  angles  of 
Interest  and  various  normalized  times.  Although  the  incident  angle,  time, 
and  roughness  condition  alter  the  magnitude  of  the  structure  function  at  the 
time  delay  t/T^  ■ 1,  the  functional  form  of  the  structure  function  for  values 
of  T/Tp<1.0  Is  unaltered  by  these  factors  even  though  the  time-resolved  LRCS 
Is  changed  significantly  (see  Section  5.2).  Clearly,  a linear  functional 
form  approximates  the  structure  function  for  values  of  t/T  <1.0  for  all  of 
the  parameters  Investigated.  Therefore,  according  to  Equations  (10)  - (12) 
the  power  spectrum  of  fluctuations  Is  of  the  form 
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Figure  13  Pulie  structure  functions  for  phenolic  carbon  cone  (6  - 0®) 
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Figure  14  Pulie  itructure  functions  for  phenolic  carbon  cone  ((-)  - 45°) 
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5. A Normalized  Cone-Structure  Functions 

Normalized  cone-structure  functions  were  calculated  using  the  Equation 
(37)  definition  for  various  incident  angles  and  roughness  conditions. 
Representative  results  for  the  more  interesting  cases  (i.e.,  large  incident 
angles)  are  shown  in  Figures  18  and  19  to  illustrate  the  effectiveness  of  the 
normalization  adopted  in  the  definition  of  the  normalized  structure  function. 
Although  the  approximate  linear  functional  form  (for  t/Tp  < 1)  noted  in  the 
previous  section  is  unaltered  by  the  normalization,  the  slopes  of  the 
structure  functions  in  general  have  been  approximately  equated  by  the 
normalization  procedure.  Thus,  it  may  reasonably  be  concluded  that  homo- 
geneous conical  targets  possess  normalized  statistical  local  homogeneity  (as 
defined  in  Section  A. 2).  Subsequent  calculations  of  intra-pulse  speckle 
effects  on  the  tlme/speckle-spectrum  signatures  of  cones  are  considerably 
simplified  by  this  fact. 
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Figure  18  Normilized  pulse  structure  functions  for  phenolic  cerbon  cone  (6  - 89°) 
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